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We present a new study of D semileptonic decays on the lattice which employs the Highly Improved 
Staggered Quark (HISQ) action for both the charm and the light valence quarks. We work with 
MILC unquenched Nf =2 + 1 lattices and determine the scalar form factor fo(q 2 ) for D — > K,lv 
semileptonic decays. The form factor is obtained from a scalar current matrix element that does 
not require any operator matching. We develop a new approach to carrying out chiral/continuum 
extrapolations of /o(q 2 )- The method uses the kinematic "z" variable instead of q 2 or the kaon energy 
E K and is applicable over the entire physical q 2 range. We find fo~* K (0) = f+^ K (0) = 0.747(19) 
in the chiral plus continuum limit and hereby improve the theory error on this quantity by a factor 
of ~4 compared to previous lattice determinations. Combining the new theory result with recent 
experimental measurements of the product /+ ^ K (0) * \V CS \ from BaBar and CLEO-c leads to a very 
precise direct determination of the CKM matrix element \V ca \, \V C3 \ = 0.961(11)(24), where the first 
error comes from experiment and the second is the lattice QCD theory error. We calculate the ratio 
f+^ K (0)/f Ds and find 2.986 ± 0.087 GeV -1 and show that this agrees with experiment. 



PACS numbers: 12.38. Gc, 13.20.Fc, 13. 20. He 

I. INTRODUCTION 

Independent determinations of each of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements and checks 
of three generation unitarity provide stringent consis- 
tency tests of the Standard Model and have become an 
important part of flavor physics. For instance, first row 
unitarity has now been checked to very high accuracy 
with \V ud \ 2 + \V US \ 2 + \V ub \ 2 - 1 = -0.0001(6) |H. Such 
precision became possible when both experiment and lat- 
tice QCD theory inputs to the determination of \V US \ 
reached subpercent levels of accuracy (the contribution 
from \V u b\ to the unitarity sum is negligible and \V u d\ is 
known with ~ 0.02% errors). 

In contrast to the elements of the 1st row, direct deter- 
minations of the 2nd row matrix elements are still much 
less precise. The latest PDG summary Q quotes 4.8%, 
3.5% and 3.2% errors for \ V c d\, \V CS \ and \ V c b \ respectively 
when one considers all ways of extracting the matrix el- 
ements. If one focuses just on determinations of \V cc i\ 
and \V CS \ from semileptonic D — > ir, Iv and D — > K,lv 
decays then until recently the total error has been at 
the 10% level and dominated by lattice QCD theory er- 
rors. Furthermore tests of 2nd row unitarity stands at 
\V cd \ 2 + \V CS \ 2 + \V cb \ 2 = 1.101 ± 0.074 and for the 2nd 
column at \V US \ 2 + \V CS \ 2 + \V ts \ 2 = 1.099 ± 0.074 0. 
In both cases the error is dominated by the uncertainty 
in \V CS \. Clearly reducing the errors in \V CS \ will have 
immediate and significant impact on flavor physics and 
on CKM unitarity tests. On the experimental front Belle 
3 , BaBar Q and CLEO-c [H| have all recently published 
precise measurements of the combination f+~* K (0)*\V CS \ 
with 3.3%, 1.4% and 1.1% errors respectively. These pre- 



cise measurements can be turned into accurate |V CS | de- 
terminations if and only if theory can provide the form 
factor fP~* K (0) with comparable precision. The main 
goal of the current work was to improve lattice QCD cal- 
culations of f+^ K (0) and the outcome is that we have 
now succeeded in reducing the theory errors from 10% 
down to 2.5%. Innovations that made this dramatic im- 
provement in errors possible include the employment of 
a new and better action for charm quarks, the use of 
an absolutely normalized hadronic matrix clement that 
does not require any operator matching, improved anal- 
ysis tools for lattice data and a new method for carrying 
out chiral/continuum extrapolations. 

One reason for the large disparity in the size of lattice 
QCD errors between kaon and D meson systems in the 
past has been the challenge of simulating quarks with 
masses as large as that of the charm quark. If "a" is the 
lattice spacing and hence ~ 1/a the cutoff in the theory, 
then for charm quarks it was thought to be difficult to 
satisfy am c = m c /cutoff <C 1. In the past this problem 
was circumvented by employing effective theories to han- 
dle the charm quark, nonrelativistic (NR)QCD, HQET 
or the "heavy clover" action of the Fermilab lattice col- 
laboration. The first Nf = 2 + 1 unquenched studies 
of D semileptonic decays on the lattice were carried out 
by the Fermilab Lattice and MILC collaborations using 
an effective theory for charm Q. This pioneering work 
predicted the shape of the form factors as a function of 
q 2 prior to subsequent verification by experiment. The 
theory errors, however, were quite large at ~ 10%. The 
calculations by the Fermilab Lattice and MILC collabo- 
rations are being improved upon and their theory errors 
should be reduced significantly soon [tJ. On the other 
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hand, working with effective theories typically leads to 
larger statistical and systematic errors than when em- 
ploying "relativistic" quark actions as can be done for 
kaon physics. This includes uncertainties coming from 
matching of heavy-light currents responsible for heavy 
meson leptonic or semileptonic decays, and from the tun- 
ing of heavy quark masses, all procedures that are much 
more complicated in effective theories. 

In recent years there has been a major shift in lat- 
tice simulations with charm quarks. With the advent of 
highly improved lattice quark actions the concerns de- 
scribed above of m c /cutoff being too large have been 
overcome and several collaborations are now working on 
charm physics without resorting to effective theories. In 
2007 the HPQCD collaboration introduced the "Highly 
Improved Staggered Quark" (fflSQ) action Q. Many 
lattice artifacts, including all 0((am c ) 2 ) discretization 
effects and all 0(a s (am c ) 2 ) and 0((am c ) 4 ) effects at 
leading order in the charm quark velocity - have been 
removed. It then becomes feasible to simulate charm 
quarks on the lattice in a fully relativistic setting with- 
out introducing large discretization errors as long as one 
works with lattice spacings a <~ 0.15fm. This last con- 
dition is easily met nowadays in typical lattice simula- 
tions. Another approach to charmed meson leptonic and 
semileptonic decays using relativistic charm quarks is be- 
ing pursued by the European Twisted Mass (ETM) col- 
laboration Q . They employ a special version of Wilson 
type quark action called the "Twisted Mass" quark ac- 
tion where discretization errors come in at 0(a 2 ). 

The HPQCD collaboration has successfully applied 
HISQ charm quarks to studies of D and D s meson lep- 
tonic decays [lOj. This formalism significantly reduced 
lattice errors in decay constant calculations. We have 
now also initiated D meson semileptonic decay studies 
based on HISQ charm and light valence quarks. In ref- 
erence [ll[ we developed and tested our approach on 
a simpler test case of a fictitious D s — > r\ Sl lv semilep- 
tonic decay (rj s refers to a pseudoscalar s — s bound 
state). In the current paper we present the first ap- 
plication of HISQ charm quarks to realistic D meson 
semileptonic decays. More specifically we calculate the 
scalar form factor fo(q 2 ) for D — > K, Iv decays. As men- 
tioned above already experiment provides us with the 
product /+(0) * IV^I- Then using the kinematic relation 
/o(0) = /+(0), an accurate calculation of fo{q 2 ) from 
the lattice leads to a precise determination of the CKM 
matrix element |V^ S |. 

We have carried out simulations on three of the MILC 
"coarse" ensembles with lattice spacing a ~ 0.12fm and 
two of the "fine" ensembles with a ~ 0.09fm. Details of 
these ensembles are listed in Table |U The five ensembles 
provide enough variation and information to allow sensi- 
ble chiral and continuum extrapolations to the physical 
world. 

In the next section we review the formalism for ex- 



TABLE I: Details of MILC configurations employed in this 
article. Ntsrc is the number of time sources used per configu- 
ration. All sea quark masses are given in the MILC collabo- 
ration normalization convention with u^ laq — (plaquette) 1 ^. 
Values for the scale variable n in lattice units, r\ /a, are taken 
from Ref. Errors in this ratio are at the ~0.1% level. 



Set 


n/a 


aUolUsea 


a Q 


^ conf 


Ntsrc 


L 3 


x N t 


CI 


2.647 


0.005/0.050 


0.8678 


600 


2 


24 3 


x 64 


C2 


2.618 


0.010/0.050 


0.8677 


600 


2 


20 3 


x 64 


C3 


2.644 


0.020/0.050 


0.8688 


600 


2 


20 3 


x 64 


Fl 
F2 


3.699 
3.712 


0.0062/0.031 
0.0124/0.031 


0.8782 
0.8788 


600 
600 


4 
4 


28 3 
28 3 


x 96 
x 96 



tracting semileptonic decay form factors from hadronic 
matrix elements of appropriate heavy-light currents. We 
describe the advantages of working with the same rela- 
tivistic action for both the heavy and the light quarks 
and explain why the form factor at q 2 = is most ac- 
curately extracted from hadronic matrix elements of the 
scalar current as opposed to from the vector current. In 
section III we introduce the HISQ action and describe 
how action parameters such as bare quark masses were 
tuned. Section IV provides further details of our sim- 
ulations including benefits derived from using "random 
wall" sources, in particular when simulating kaons with 
nonzero momenta. Section V describes our fitting strat- 
egy. We have invested considerable effort into developing 
improved fitting methods in order to extract the form fac- 
tors of interest with subpercent errors for each ensemble 
and for all kaon momenta needed to cover the physical 
q 2 range. 

In section VI we take the form factor results from the 
five ensembles and extrapolate to the chiral/continuum 
limit. To this end we have developed a new extrapola- 
tion method that can be used over the entire physical q 2 
range, (Md — Mk) 2 > q 2 > 0. Since we are interested 
in the form factor at q 2 = 0, it is important that any 
extrapolation scheme work all the way down to q 2 = 
where the kaon in the D meson rest frame has energy 
Ek ~ lGeV. The new approach uses the "z-cxpansion" 
of Refs. [l3-[l5| to parameterize the kinematics (the de- 
pendence on Ek or cquivalently on q 2 ). The coefficients 
of this expansion are then allowed to be functions of the 
light and strange quark masses and of the lattice spac- 
ing. We find that good fits to all our data are possible 
with such an ansatz and the resulting /o(0) in the chi- 
ral/continuum limit is very stable against higher order 
terms in this ansatz. 

Section VII summarises our final results for /o(0) = 
/+(0) at the physical point and explains our error bud- 
get. We incorporate experimental input from BaBar Q 
and CLEO-c [f| to determine the CKM matrix element 
| V cs | and compare with values listed in the PDG and with 
expectations from CKM unitarity. In Section VIII we 
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collect results obtained as "side products" of our anal- 
ysis of semileptonic decay three-point hadronic matrix 
elements, namely results coming from two-point correla- 
tors such as decay constants of the pion, kaon, D and 
D s mesons. These two-point results provide nontrivial 
tests of our mass tunings, fitting and chiral/continuum 
extrapolation strategies leading to greater confidence in 
our form factor (i.e. three-point) results as well. Finally 
section IX gives a summary and addresses future plans. 
Several appendices cover details of Bayesian fits employed 
in this article. And in Appendix C we carry out the chi- 
ral/continuum extrapolation of /o(<? 2 ) using an approach 
that differs completely from the z-expansion method of 
section VI and is based instead on Chiral Perturbation 
theory. We show that the two extrapolation methods 
give results in very good agreement with each other. 



II. FORMALISM 

To study the process D — » K, Iv one needs to evalu- 
ate the matrix element of the charged electroweak cur- 
rent between the D and the K meson states, (K\{V^ — 
A^)\D). Only the vector current U M contributes to the 
pseudoscalar-to-pseudoscalar amplitude and the matrix 
element can be written in terms of two form factors 
f+{q 2 ) and fo(q 2 ), where = — p^ K is the four- 
momentum of the emitted W-boson. 



TABLE II: Action parameters and rj c and rj a masses. 1 + e 
enters the coefficient of the Naik term for charm propagators. 



Set 


araoc 


amos 


amoi 


1 + e 


aM 7)c 


aM Vs 


CI 
C2 
C3 


0.6207 
0.6300 
0.6235 


0.0489 
0.0492 
0.0491 


0.0070 
0.0123 
0.0246 


0.780 
0.774 
0.778 


1.7887(1) 
1.8085(1) 
1.7907(1) 


0.4111(2) 
0.4143(2) 
0.4118(2) 


Fl 
F2 


0.4130 
0.4120 


0.0337 
0.0336 


0.00674 
0.01350 


0.893 
0.894 


1.2807(1) 
1.2751(1) 


0.2942(1) 
0.2931(2) 



We have checked the feasibility of this renormalization 
scheme and extracted preliminary Z t and Z s values for 
the test case of D s — > -q si Iv mentioned above in Ref . [TlT] . 
We plan to apply this fully nonperturbative renormaliza- 
tion scheme to evaluate (it\V^\D) and (K\V^\D) relevant 
for realistic D — > n, Iv and D — > K, Iv semileptonic de- 
cays in the near future. In the present article, however, 
we will focus on the form factor f+(q 2 ) just at q 2 = 0, 
since this is all that is needed to extract |V CS |. We do 
this by exploiting the kinematic identity /+(0) = /o(0), 
and concentrating on determining the scalar form factor 
fo(q 2 ) as accurately as possible. The best way to proceed 
is to evaluate the hadronic matrix element of the scalar 
current rather than of the vector current. From eq.@ 
one then has, 



(K\V\D) 



= f?^ K (q 2 ) 



M 2 



■tic — 



M 2 D 



Ml 



(1) 



with = g-f^^c. As described below, we find it useful 
to consider also the matrix element of the scalar current 
S = 



(K\S\D) = 



M 2 - M 2 
iu D iu K D _ 

Jo 

moc - m 0s 



■K 



(<7 2 )- 



(2) 



In continuum QCD one has the PCVC (partially con- 
served vector current) relation and the vector and scalar 
currents obey, 



(m 0c ~m 0s )(S cont -). (3) 

In fact PCVC is the reason why the same form factor 
fo > ^ K (q 2 ) appears in eqs.([T]) and ([2]). On the lattice it is 
often much more convenient to simulate with vector cur- 
rents vE'qi7 ai 4'q2 that are not exactly conserved at finite 
lattice spacings even for Ql = Q2. Such non-exactly- 
conserved currents need to be renormalized and acquire 
Z-factors. We are able to carry out fully nonperturbative 
renormalization of the lattice vector current by imposing 
PCVC. In the D meson rest frame the condition becomes, 



(M D -E K )(V tatt -)Z t +p K -(V latt )Z s = (m 



-m 0s )(S latt -) 
(4) 



(m 0c - m 0s ){K\S\D) 



M 2 D 



M 2 



(5) 



The numerator on the right-hand-side is a renormaliza- 
tion group invariant combination. This is true even in our 
lattice formulation, because we use the same rclativistic 
action for both the heavy and the light valence quarks. 
Moreover, cq. ([5]) allows a lattice determination of fo(q 2 ) 
and hence also of /+(0) = /o(0) without any need for 
operator matching. Using eq. (|5|) and going to the con- 
tinuum limit is straightforward, because our action is so 
highly improved even for heavy quarks. 



III. THE HISQ ACTION AND TUNING OF 
ACTION PARAMETERS 



The HISQ action was introduced in Ref.[8| and repre- 
sents the next level of improvement of staggered quarks 
beyond the AsqTad [ll[ action. Relative to the latter, 
the HISQ action reduces taste breaking effects by ap- 
proximately a factor of three on MILC coarse and fine 
lattices. For heavy quarks such as charm the HISQ ac- 
tion includes one adjustable parameter e which modifies 
the "Naik" term already present in the AsqTad action 
|a 2 A^ ii^a 2 A^. In Ref.0 e was adjusted nonper- 
turbativcly to get the correct dispersion relation for the 
r] c meson. It was found that one ends up with e's close to 
estimates coming from requiring that the tree-level quark 
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propagator have speed of light c(p) = 1. In the current 
simulations we have set e in the charm propagators equal 
to its tree-level value. The size of e decreases rapidly with 
decreasing am qi so for strange and light quarks e can 
safely be set equal to zero. Table HIl lists the bare valence 
masses and the e(charm) values employed for the five en- 
sembles. The bare charm and strange quark masses were 
tuned in order to obtain the correct r\ c and r\ s meson 
masses, ara^i was chosen so that the ratio of HISQ va- 
lence masses rriQi(valence) /mo s (valence) ended up close 
to the analogous ratio of light and "physical" strange 
quark masses for AsqTad quarks on the same ensemble 1 . 
This makes the valence light quark mass approximately 
equal to the sea light quark mass, both measured rela- 
tive to a physical, tuned strange quark mass. Figures [1] 
and [5] show our results for the r\ c and r\ s masses on the 
five ensembles together with the target values. One sees 
that tuning has been achieved very accurately. For r\ c 
our target value is M*^ e * = 2.9852(34)GeV which dif- 
fers slightly from the experimental M^ x J er - = 2.9803GeV 
since we adjust for the absence of electromagnetic and 
annihilation effects in the lattice simulations The 
target value for the rj s is M Vs = 0.6858(40)GeV [l9|. 
The data points in Figs. [T] and [2] show only statistical 
and ri/a errors. We first determine r\ x M meson using 
the precisely know r\ja for each ensemble. At this point 
our data points have errors coming from both statistics 
and from the ~0.1% uncertainty in r\/a (with the latter 
dominating). For the purposes of using a physical scale 
on the vertical axis we then convert r% x M meson to MeV 
using rl = 0.3133(23)fm (rf 1 = 0.6297(46)GeV) Ql]. 
For reasons explained below, we find it more informative 
not to include the ~ 0.7% uncertainty in the physical 
value of r\ in these plots. Including this error will affect 
all five data points in the same way without changing 
relative uncertainties. 

Once the bare charm and strange quark masses have 
been fixed for each ensemble there are no adjustable pa- 
rameters left when one goes on to determining other me- 
son masses such as Mu or Mjj s , decay constants fo, 
Id,, fx etc. or scmilcptonic form factors. In Fig. [3] 
we show results for Mjj and Md s - Again the errors on 
the data points reflect statistical and r\/a errors only. 
Any changes in the physical value of r\ will shift all 
five data points uniformly without affecting their relative 
positions. Furthermore chiral/continuum extrapolations 
would be carried out on r\ x M meson with r\ and its error 
coming in only after having extracted the physical limit. 
By omitting the full r\ errors in Fig. [3] one can more eas- 
ily identify discretization effects and light quark mass de- 
pendence. For instance, for Md s the difference between 
the coarse and fine ensemble results is at the ~ 6MeV 
level or ~ 0.3% and the sea light quark mass dependence 
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FIG. 1: Tuning of the charm quark mass via the r\ c meson 
mass. Errors on the simulation results include statistical plus 
the ~ 0.1% errors coming from r\/a. The smaller black error 
bars on the lattice data indicate the statistical errors. The 
"experimental" r\ c mass has been adjusted to take into ac- 
count the lack of annihilation and electromagnetic effects in 
our lattice calculation. 



is essentially nonexistent. The 0.3% discretization effect 
should be compared to the ~ 0.7% uncertainty in r\. One 
lesson to be learnt from this is the importance of tuning 
quark masses accurately enough so that results on the 
different ensembles agree to within the smaller r\ja er- 
rors and not just to within the larger n error. Otherwise 
it would not be possible to have data points lying along 
smooth curves as in Fig. [3] where discretization and light 
quark mass dependence can be clearly identified and dis- 
tinguished from mass tuning and r± errors. Based on 
Figs. [U [2] and [3] we believe the HISQ action parameters 
have been fixed accurately enough in preparation for go- 
ing on to D semileptonic decays. Further consistency 
checks, such as determinations of decay constants will be 
given in section VIII. 



IV. SIMULATION DETAILS 

The goal is to determine the hadronic matrix element 
(K\S\D) in eq.((5]) via numerical simulations. The start- 
ing point is the three-point correlator, 



3 ipK-(z—x) 



x y 



($ K (x,t )S(z,t) & D (y,t -T)). 



(6) 



$|j and 3>i<- are interpolating operators that create a D 
meson or annihilate a kaon respectively and S = a 3 S is 
the scalar current in lattice units. We also work with 
dimensionless fermion fields. Eq. © corresponds to first 



1 Note that HISQ and AsqTad masses are not numerically the 
same, since the two actions arc different. 



5 



700 



690 



680 



670, 



0.6% error 



I © $ 1 



□ coarse lattice 
fine lattice 
O Target 



0.1 0.2 0.3 0.4 0.5 

m,/m 



FIG. 2: Tuning of the strange quark mass via the rj a 



for each configuration we obtained results for several T 
values, T = 15 and 16 on coarse and T = 19, 20, and 
23 on fine lattices. We will see later that having data at 
many T values significantly reduces errors in extracted 
three-point amplitudes. To further improve statistics, for 
each to and T value we also evaluated the time-reversed 
three-point correlator, essentially the same as cq.© but 
with acting on time slice to+T and the scalar current 
inserted at T + to > t > to- 

As is well known, with the HISQ action, each flavor 
of quark comes in N tas te copies called "tastes" . One has 
Ntaste = 4 when working with one-component staggered 
fields and Ntaste — 16 for four-component "naive" fields. 
In the naive fields language the interpolating operators 
an d the scalar current S become, 



1- 



D - 4*c75*/, 



(7) 



and 



5" = * s * c 



(8) 
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FIG. 3: The D and D s meson masses. As explained in the 
text although the vertical axis uses physical units, the data 
points show statistical and n /a errors only. The experimental 
points have not been adjusted for electromagnetic effects that 
are absent in the lattice simulations. Such effects have been 
estimated to be less than 0.2% 12011. 



creating a zero momentum D meson at time to~T which 
propagates to time slice t with to > t > to — T . At 
time slice t the scalar current S = ^ s ^c converts the 
charm quark inside the D meson into a strange quark 
while inserting momentum pk- The resulting kaon is 
then annihilated at time to- 

In our simulations we have picked to the location of 
the kaon operator <&k randomly and differently for each 
configuration in order to reduce autocorrelations. We 
also worked with N tsrc values of to per configuration (see 
TablclJ) placed N t /N tsrc time slices apart (N t is the total 
number of time slices for our lattices). Once to was fixed, 



These are all single site bilinears. ^d/k correspond to 
taste non-singlet "Goldstone" pseudoscalars and the fac- 



tors of \ 



serve to divide out traces over taste 



space. The scalar current in eq.© is a taste singlet 
current. Carrying out the contractions over fermionic 
fields in C 3pnt and using the well known relation between 
naive quark propagators G^(x,y) and one-component 
field propagators G x (x,y), 



with 



G*(x,y) = n(x)rt(y)G x (x,y) 



fi=0 



(9) 



(10) 



one obtains, 



(^ K (x)s(z)^ D (y)) 

—Tr {Gq,, s (x, z)G*. c (z, y)^Gin^{y, x)j 5 } 
16 

1 

16 



xG XtS (x, z)G XyC (z, y)G Xi i(y, x)} 



= -£<t>(y)4>(x)tr {G x , s (x, z)G XtC (z, y)G X; i{y, x)} 

= -^{v)4>{.z)tr {G XtS (z, x)G XtC (z, y)G Xi i(y, x)} . 

(11) 

V y 

4>{y) stands for (—1)^^ and similarly for cj>(x) and 
4>{z). In the last step we have used G x (x,y) = 4>{x) * 
<f>(y)Gt(y,x). "Tr" is the trace over spin and color and 
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"tr" the trace only over color. The three-point correlator 
can now be written as 



c^ nt (t , t, t, p K ) = -1 e E E ^ k ' (S ~ s 



x y z 



-cf>(y)4>(z){tr {G\ ltS (z,x)G XtC (z,y)G Xt i(y,x)}), 

(12) 

with xq = to, yo = to — T and zq = t and where () now 
stands for average over configurations. One sees from 
cq.(fT2"j) that strange HISQ propagators are needed going 
from (x, to) to general z and light propagators again from 
(x, to) to general y. If one actually wanted to carry out 
the 7j one would need a strange and a light propa- 
gator from each spatial point on time slice to and that 
would be prohibitively expensive. A common approach 
is to give up on doing the J2s an< ^ to use "local sources" 
where x is fixed at some xo, e.g. x = 0. One then has, 

c?:: t (t ,t,T,p K ) = j2J2 eti!K '* x 
i 



V z 



'-6(y)(j)(z)(tr {G x (z,xi 0C )G x>c (z,y)G Xi i(y,xi 0C )}), 



(13) 



with xioc = (0,to). Momentum conservation will ensure 
that only kaons with momentum px contribute and be 
picked out at time slice to- On the other hand "random 
wall" sources allow us to carry out the jj ^ without 
having to invert at each spatial point. This can be seen 
by writing, 



C 3 Z nt (t ,t,T, 



PK) = 



73EEEE 1 



,ipK-(z—x) 



x x y 



{ ct>{y)<t>{z){tr{G X!S {z, x)G x Jz, y)G xJ (y, x') 



(14) 



t;(x) is a field of random U(l) phases, and (£* (x)^(x' )) = 
Sg^/ ensures that (fT4")l reduces to (|T2"j) after averaging over 
gauge field configurations. So in the random wall source 
approach one calculates the light quark propagator G Xi z 
using the source y=f an d G XyS , the strange 

propagator by using the source ^= Ej?P) e * x - Only 
one light quark inversion is required in addition to a sep- 
arate strange quark inversion for each px- In this way 
one obtains the random wall propagators, 



G r x %y,t ) 



and 

G^ s (z,t ;p K ) 



1 



V T-' 



EG^(2/,*W) 



EG X)S (z,z)£(af)e* 



(15) 



(16) 
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FIG. 4: Comparison between local and random wall sources 
for p = £(0,0,0) (left plot), and for p = £(1,1,1) (right 
plot). 



The expression for the three-point correlator becomes, 

c 3 r r\t , t, r,p K ) = E E e%pK ' s x 
1 



V z 



^( y )^z)(tr{G x w J( Z ,to;p K )G x Az,y)G x w l ( y ,to)}), 

(17) 

The charm propagator in ea. (|17p is obtained by inverting 
at time slice yo = t — T with source J2y ( l ) (y)G x w i(y,to)- 
In this way one gets the "sequential" charm propagator, 

G^(z, t , T) = E <KV)G X A*> v)G x w l {y, t ), (18) 



and 



C 3 Z nt (t ,t,T, 



PK) 



WK-Z x 



-<t>(z)(tr{G x w Xz,t \PK)G x %z,to,T)}). (19) 

The most costly part of our simulations is calculating 
the random wall strange quark propagators of eq. (fTB")) . A 
separate inversion is required for each px (e.g. 8 inver- 
sions for the different combinations (±1,±1,±1)). On 
the other hand, when we change the T values only G x f 
of eq. (fl8|) needs to be recalculated and one inversion suf- 
fices for all momenta. This is one of the reasons why the 
full kaon momentum px is put into (|16p and none into 



In Fig. U we show comparisons of percentage errors in 
three-point correlator data of local sources versus random 
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wall sources. One sees significant improvement coming 
from random wall sources. These tests were carried out 
in the test case D s — > rj Sl lv calculations and with less 
than the full statistics. For D — > K,lv we immediately 
went to random wall sources. 

In addition to the three-point correlators, as we will 
show in the next section, several two-point correlators are 
needed in order to extract the matrix element (K\S\D). 
They are, 

c 2 D mt (t, to ) = ±j2T,(*°(y^D&to)), (20) 



and 



1 EE e 

v 



L 3 



^-y\^ K ^t)^ K {xM))- 



(21) 



The -p- ^ - can again be implemented via random wall 
sources. As already mentioned in the previous section we 
also calculated correlators for the r/ c , r) s and D s mesons in 
exactly the same way as C^J 1 ^- in order to carry out and 
check mass tunings. We have accumulated simulation 
data for the three- and two-point correlators described 
in this section for the ensembles of Table HI In the next 
section we explain how hadronic matrix elements such 
as (K\S\D) and meson masses and decay constants are 
extracted from this data. 



V. FITS AND DATA ANALYSIS 

The interpolating operators $d and $k do not cre- 
ate just the ground state D meson or kaon that we are 
interested in but also excited states with the same quan- 
tum numbers. With staggered quarks there is the further 
complication that in addition to regular states so-called 
"parity partner" states can contribute, whose energies are 
measured relative to iir so that e~ Et — > (— l) t e~ Et . The 
D meson correlator, for instance, has the t dependence 
(we set to = for simplicity), 

N D -1 

C 2 D pnt (t) = ]T ^(e-^+e-*?^-*)) 

N'n-l 

+ ]T ^(-l)*( e -^t + e -^(%- t)) _ 



fe=0 



(22) 



We are interested in the ground state D meson contribu- 
tion with amplitude, 



l n_\{*D\D)\ i 



2M D a 3 



(23) 



Similar relations apply for other mesons. Only in the case 
of equal mass mesons (n, i] s or rj c ) at zero momentum 



are the oscillatory contributions absent. The three-point 
correlators such as ea. (|19|) will have contributions from 
regular and oscillatory states for both the kaon and the 
D meson. The rather complicated t and T dependence 
is then given by, 

N K -i n d -i 

C 3 P"'(i,T)= £ A jk e- E ? t e- E ?<?- 1 > 

j k 

N K -lN' n -l 

+ E E Bjke-^e-^^i-l)^ 

3 k 
N' K -lN D -l 

+ E E C jk e- E ^e-^^(-iy 

j k 
N' K -lN' n -l 

+ E E D jk e- B '* t e- B '* DlT - t \-l) t (-l) lT - t l 

j k 

(24) 

We will only consider the region < t < T and take 
T « N t so that any contributions from mesons prop- 
agating "around the lattice" due to periodic boundary 
conditions in time can be neglected. The relevant ampli- 
tude here is, 



An 



($k\K) (K\S\D) (D\$ d ) 
(2E K a 3 ) (2M D a 3 ) 



(25) 



From cqs. ([2"3"|) and ([S]) one sees that our sought after 
hadronic matrix element is given by, 



(K\S\D) = 2y/M D E K 



A 



(in 



(26) 



So our goal is to extract the combination on the right- 
hand-side of (|26[) as accurately as possible and with any 
correlations among the errors of the individual compo- 
nents, Ann, &o °, Md and Ek taken properly into ac- 
count. 

We have carried out simultaneous fits to C^ nt , C 2 ^ nt 
and the three-point correlators with different T values, 
C 3pnt (t,Ti) i = 1,2.... following the fit ansaetze of 
eqs. (|2"2"|) and (|2"4")l . Two (three) different T values are used 
for the coarse (fine) ensembles. This allows us to evalu- 
ate (|2l)|) within one fit. The two-point correlators were fit 
for lvalues between i TO j„ = 2(2) and t max = 30(20 ~ 30) 
for the K and D respectively for coarse lattices. For fine 
lattices, lvalues were used between t m i n = 2 ~ 4(2) and 
tmax — 30(30). For the three-point correlators we used all 
the data between t = 2 and t = T — 2 for coarse lattices, 
and t = 3 and t = T—2 for fine lattices. Simultaneous fits 
with multiple T and taking different fit ranges for differ- 
ent correlators were also helpful to reduce the statistical 
errors, since it allows us to extract maximum information 
from both the three-point and two-point correlators. The 
number of exponentials in our fit ansatz was varied to test 
for stability of fit results. For our final fits we ended up 
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N K {Nd) is fixed at 3 (4). 
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FIG. 6: aMn versus Nd- Green circles are from simultaneous 

C 2pnt _ C 3pnt fitg ^ diamond jg from fitg tQ j ugt the 



choosing around 3^4 for No and Nk- N' d , k was taken 
to be mostly N D / K — 1. Figs [5] - [9] show some results for 
A o, Md, bP, Ek, bff versus Nd or Nk for ensemble CI 
at kaon momentum p = (0, 0, 0). All our fits are carried 
out using Bayesian methods [21( • We describe choices for 
priors and prior widths in the Appendix. 

We have found that using data from several C 3pnt 
with different T-valucs helps greatly in reducing statis- 
tical/fitting errors. Fig. [10] compares results for fo(q 2 ) 
for ensemble C2. One sees that having two rather than 
just one C 3pnt (t,Ti) involved in the simultaneous fit re- 
duces errors and that this effect is most pronounced when 
one combines an even T with an odd T. It may not be 
surprising that improvements are achieved from multi-T 



fits. From the fit ansatz ([24]) one sees that having more 
C 3pnt, s does not increase the number of fit parameters 
(Ajk etc.) although the amount of data and hence of 
information given to the minimizer is increased. Fur- 
thermore if Ti + T 2 is odd, then (-l) Tl ~* and (-l) 7 " 2 " 1 
have opposite signs and C 3pnt (t, Ti) and C 3pnt (t, T 2 ) will 
provide more independent information. 

Most of the fit parameters such as Md or b$ etc. that 
one gets from the simultaneous C 2pnt - C 3pnt fits can 
also be determined by just fitting the two-point correla- 
tors by themselves. Results from the C 2pnt fits are also 
shown on Figs. [6] -[9] and provide consistency checks. One 
interesting outcome is that two-point correlator param- 
eters are more accurately determined via simultaneous 
fits with three-point correlators than when they are fit 
alone. This is especially noticeable for D meson corre- 
lators, namely for Md and bP and has implications for 
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determinations of the decay constant Jd- 
related to through, 



afi 



m . 



m i 



M D 



l2bl_ 
aM D 



The latter is 



(27) 



In Figs. [TT] and [T^] we compare results for Mr> and fjj 
using cither pure two-point fits or simultaneous fits. One 
sees the significant improvement coming from the simul- 
taneous fit. Doing simultaneous fits gives a better handle 
on excited state contributions because they contribute 
differently to two-point and three-point correlators. This 
is effectively similar to adding smearings to the correla- 
tors. In section VIII we will discuss extracting fjj in the 
chiral/continuum limit. Having simultaneous fit results 
will make this determination more accurate than what 
can be achieved from pure two-point correlators. This 



FIG. 12: 

fits 



Comparison of fo from two-point and simultaneous 



appears to be a bonus side product of semileptonic decay 
studies. 



In Tables IIII1 IIVI and [V] wc summarize our main fit 
results. One sees from Table [V] that we were able to 
determine /o(pW) with errors ranging from ~ 0.2% at 
zero momentum to ~ 0.9% at our highest momentum. 
In Fig. [13] we plot the square of the "speed of light" c 2 (p) 
for the kaon. 



E%{p)-M% 



P 



(28) 



One sees that the relativistic dispersion relation is satis- 
fied to about 1 ~ 2%. Wc will check the effect of devi- 
ations from exact continuum dispersion relations on our 
final results for form factors in later sections. 
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TABLE III: Fit results for aM D , af D , aM Ds and af Ds . Num- 
bers in brackets are from two-point fits, whereas the rest come 
from simultaneous fits. 



TABLE V: Fit results for f (p f 



Set 


aM D 




aM Ds 


a fo s 


CI 


1.1393(7) 


0.1372(4) 








(1.1398(21)) 


(0.1373(18)) 


(1.1876(5)) 


(0.1539(6)) 


C2 


1.1595(8) 


0.1423(4) 








(1.1574(19)) 


(0.1405(17)) 


(1.2008(6)) 


(0.1560(7)) 


C3 


1.1618(5) 


0.1464(3) 








(1.1620(10)) 


(0.1464(9)) 


(1.1899(5)) 


(0.1553(4)) 


Fl 


0.8141(3) 


0.0971(2) 








(0.8152(8)) 


(0.0979(5)) 


(0.8473(2)) 


(0.1083(2)) 


F2 


0.8197(3) 


0.1007(2) 








(0.8191(5)) 


(0.1005(4)) 


(0.8435(2)) 


(0.1078(2)) 



TABLE IV: Fit results for clMk, o,Ek(p), aftc, am, and afiv- 
Numbers in brackets are from two-point fits, whereas the rest 
come from simultaneous fits. 



Set 


clMk 


aE K 


aE K 


aE K 






(1,0,0) 


(1,1,0) 


(1,1,1) 


CI 


0.3122(2) 


0.4081(7) 


0.4837(9) 


0.5469(20) 




(0.3122(2)) 


(0.4083(7)) 


(0.4842(9)) 


(0.5477(24)) 


C2 


0.3285(5) 


0.4531(16) 


0.5525(17) 


0.6373(32) 




(0.3285(3)) 


(0.4536(13)) 


(0.5532(18)) 


(0.6382(33)) 


C3 


0.3572(2) 


0.4750(9) 


0.5720(10) 


0.6524(22) 




(0.3572(2)) 


(0.4755(9)) 


(0.5722(11)) 


(0.6524(35)) 


Fl 


0.2285(2) 


0.3203(7) 


0.3919(9) 


0.4559(15) 




(0.2286(2)) 


(0.3185(12)) 


(0.3896(23)) 


(0.4506(47)) 


F2 


0.2460(1) 


0.3340(4) 


0.4014(7) 


0.4609(11) 




(0.2458(2)) 


(0.3334(7)) 


(0.4015(10)) 


(0.4616(16)) 



Set 


afi< 




afw 


CI 


0.1011(1) 








(0.1011(1)) 


(0.1599(2)) 


(0.0893(1)) 


C2 


0.1044(1) 








(0.1045(1)) 


(0.2108(2)) 


(0.0949(1)) 


C3 


0.1079(1) 








(0.1079(1)) 


(0.2931(2)) 


(0.1023(1)) 


Fl 


0.0721(1) 








(0.0721(1)) 


(0.1344(2)) 


(0.0645(1)) 


F2 


0.0748(1) 








(0.0747(1)) 


(0.1873(1)) 


(0.0697(1)) 



VI. CHIRAL AND CONTINUUM 
EXTRAPOLATIONS USING THE ^-EXPANSION 

The twenty entries in Table [V] summarize our results 
for the form factor fo(q 2 ) evaluated on the five ensembles 
of Table |I] with four different momenta pK (including zero 



Set 


fo (0,0,0) 


/o(l,0,0) 


/o(l,l,0) 


/o(l,l,l) 


CI 
C2 
C3 


1.022(3) 
1.023(4) 
1.010(2) 


0.916(3) 
0.885(5) 
0.883(3) 


0.846(3) 
0.807(4) 
0.803(4) 


0.794(6) 
0.758(7) 
0.754(5) 


Fl 
F2 


1.019(2) 
1.011(1) 


0.876(3) 
0.874(2) 


0.796(3) 
0.792(2) 


0.745(4) 
0.739(4) 



1.04 



1.02 



0.98- 



0.96 - 




FIG. 13: Speed of light squared versus momentum for the 
kaon from the 2-pt fit 



momentum) per ensemble. The kaon energy in the D rest 
frame, Ek, is related to q 2 via, 



Ml + M K - 2M D E K , 



(29) 



and the physical region is < q 2 < q 2 ,^ = (Mjj — Mk) 2 ■ 
Mk, Ek and Md are given for the different ensembles 
in Tables IIIII & IIVI The next step is to extrapolate 
the data of Table |V] to the chiral/continuum physical 
limit. As is well known, chiral extrapolations for form 
factors are much more subtle than for static quantities 
such as masses or decay constants. The main reason 
for this is that form factors depend not only on me- 
son/quark masses but also on a kinematic variable such 
as q 2 (or cquivalently on Ek)- Kinematic variables are 
themselves functions of meson masses. It is not suffi- 
cient to parameterize just the light quark mass depen- 
dence of form factors. One must at the same time cap- 
ture the kinematic variable dependence correctly for each 
value of the light quark mass and the lattice spacing 
(i.e. for each of our ensembles). Furthermore we are 
interested in a parameterization that works over the en- 
tire physical kinematic range. In the chiral limit, q 2 



ranges between < q 
0.495 < E K < l.OGeV. 



< 1.9GeV and Ek between 
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The z-expansion 

In addition to q 2 and Ek a third kinematic variable 
has proven to be convenient in semileptonic form factor 
studies, in p articular in recent analysis of B — > 7r, Iv de- 
cays [Mill- 



z{q ,t ) 



y/t+ -q 2 -VU 



to 



q 2 + y/t+ - t 



(30) 



where t± = (Md ± Mk) 2 and to is a free parameter that 
defines the zero of the z-variable, z(q 2 = to, to) = 0. By 
going to the z-variable one is mapping the cut region 
t+ < q 2 < oo in the complex q 2 plane onto the circle 
\z\ = 1 and — oo < q 2 < t + onto z G [—1, 1]. The physical 
region < q 2 < q 2 lax = t- corresponds then to an even 
smaller region around z = 0. For instance for the choice 
to = 0.5t_ and for physical values of Mo and Mk one 
has —0.057 < z < 0.046. In other words one always 
has \z\ < 0.06 in the physical region and this should 
make z a good variable for a power series expansion. As 
discussed in the literature using analyticity properties of 
form factors one can write, 
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FIG. 14: / (0) versus k max for ensembles C2 and Fl. P(q 2 ) 
(1 — q 2 /{Mo* ) 2 ) everywhere except for at k 



,x — 2 where 
2\ „_„ 



results for both P(q 2 ) = 1 and P(q 2 ) = (1 - q 2 /(M D * o 
shown. 



fo(q 2 ) = 



1 



P{q 2 )<S>o{q 2 M) 



J2a k (t )z(q 2 ,t ) k . (31) 



k=0 



The function P(q 2 ) in the denominator is there to fac- 
tor out any isolated poles in the region t- < q 2 < t+ 
below the DK threshold at q 2 = t + . In the case of 
D — >• K semileptonic decays the charm-strange scalar 
current has the same quantum numbers as the -D* (2317) 
+ meson so that one choice for P(q 2 ) would be P(q 2 ) = 
(1 - q 2 /(Al D * o ) 2 ). We have worked with both P(q 2 ) 
: 1 - q 2 /{M D * Q ) 2 ) and P{q 2 ) = 1 and find that although 
the expansion coefficients a k depend on the choice for 
P(q 2 ), in either case the data can be reproduced very 
well with just a few terms (as we will discuss below, with 
just three terms) in the z-expansion. 

For the "outer function" $0 we adopt the choice given 
inRef. [l|. 



$o0? 2 ,to) 



3M_ 



(v^ 



32ttxo 
{t + -q 2 f/ 2 {^T + 



f s/t + - t ) x 

Y + VtT^) 1/2 



(32) 



Xo fias been calculated in the literature using QCD per- 
turbation theory and the Operator Product Expansion 
(OPE). An expression including 0(a s ) and condensate 
contributions is given for instance in Ref. [l4[ ■ For fixed 
charm quark mass xo is a constant and affects just the 
overall normalization of the dfc's. For simplicity we ignore 
the condensate contributions, which are of 0(m~ 3 ) and 
0{m~ A ) respectively, and retain just the tree-level and 
0(a s ) contributions. Any other choice would just mean 



a common overall rescaling of the expansion coefficients 

Ofc. 

Testing the z-cxpansion with individual fits 

In order to test the usefulness of the z-expansion we 
first fit fo{q 2 ) separately for each individual ensemble us- 
ing the ansatz of ea. (|3T1) with J2k ^ J2k=o"- ^ n eacn 
case values used for Mb, Mk and hence also for t± 
were those specific to that ensemble. We employed a 
common choice for to, to = 0.942GeV 2 corresponding to 
to = 0.5 x £c° ntmuun \ As is well known from z-expansions 
in general, other choices for to made no difference in re- 
sulting fit curves. We find that good individual fits are 
possible once k max reaches k max = 2 and that fit curves 
are then very stable with respect to further increases in 
kmax- In Fig. [14] we show representative results for two 
ensembles C2 and Fl for fo(q 2 = 0) versus k max . We 
also show results for two different choices, P(q 2 ) = 1 and 
P(q 2 ) = (1 - q 2 /{Ml,, Q ) 2 ). One sees that fit results are 

very insensitive to these changes in k max or P(q 2 ). At 
finite lattice spacing we have used Md« = M l ^ tlce + SM 
with SM = [M D * Q (0 + ) - M Ds (0-)] exper .. These individ- 
ual z-expansion fit tests demonstrate (as advocated in 
the literature) the efficiency of z-expansions in capturing 
the kinematics of form factors with just a small number 
of parameters and in a model independent way. 

In Fig. [T5]we give examples of z-expansion fits to indi- 
vidual ensembles, specifically for ensembles C2 and Fl. 

Simultaneous modified z-cxpansion fit 

Having verified the efficacy of the z-expansion in fits 
to individual ensembles we turn next to modifying the 
fit ansatz to enable extrapolation to the physical limit. 
All kinematic properties that depend on q 2 are absorbed 
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FIG. 15: Individual z-expansion fits for ensembles C2 and 
Fl. 



TABLE VI: The expansion coefficients <Zi *Di, i = 0,1,2 from 
a simultaneous z-expansion fit to all data. 



Set 


a * D 


ai * Di 


0,2 * D2 


CI 


0.095(2) 


0.085(21) 


-0.07(11) 


C2 


0.094(2) 


0.083(21) 


-0.07(11) 


C3 


0.089(2) 


0.079(21) 


-0.07(11) 


Fl 


0.094(1) 


0.084(18) 


-0.07(11) 


F2 


0.091(1) 


0.080(17) 


-0.07(11) 


physical 


0.097(2) 


0.088(18) 


-0.07(11) 


limit 









by P, $0, and z. A natural way to distinguish between 
ensembles is to let ak — > a,k * Dk , where Dk contains the 
light quark mass and lattice spacing dependence as shown 
below (we set k rnax = 2 and P(q 2 ) = (1 - q 2 /{M D >J 2 )). 

fo(q 2 ) = 



1 



(a D + a\DiZ + a 2 D 2 z' 



P(q 2 ) $0 
x(l + b 1 (aE K ) 2 +b 2 (aE K ) i ), 



(33) 



where, 

Di = 1 



+€ 



c\xi + c 2 Sx s + c\xilog(xi) + di(am c 
iK) 4 !/, \^8Ml + 5M 



Xl 

6x s 
SM 2 



Ml 

K 



M 2 phl 



(47T/.) 2 

_}_ {{M A sqTadf _ {M HIS Q? 



(34) 
(35) 

(36) 
(37) 



FIG. 16: Chiral/continuum extrapolation of /o(<7 2 ) versus 
E\ from the modified z-expansion ansatz. The data points 
are coarse lattice points. Three individual curves and the 
extrapolated band are from a fit to all five ensembles. 



8M 2 K 



((M 



AsqTad 
K 



{M HISQ 



K 



?)■ (38) 



In cq. 34, we put typical analytic terms for light valence 
(xi and 5x s terms) and sea quark mass (SM n and SMk 
terms) dependence. We quote M^ sqTad and M^ sqTad 
from rcf. [l^- For the chiral logs, we only include 
up/down quark contributions. The strange quark chiral 
logs are close to a constant that can be absorbed into the 
Oi's. There are two distinct sources of lattice spacing de- 
pendence. (am c ) 2 and (am c ) 4 terms are due to the heavy 
quark discretization error, and (aEx) 2 and (ai?if) 4 terms 
are introduced to estimate the discretization errors due 
to finite momentum. Since we want the ciiDi to be in- 
dependent of the momentum, the ciEk terms are placed 
separately outside the ^-expansion. We include lattice 
spacing dependent terms up to fourth power, however 
we tested with even higher terms and confirmed that the 
higher terms are negligible (see Fig.l2"2"j). We have carried 
out simultaneous fits to all the data of Table IVl using the 
above ansatz and find that very good fits are possible. 
Figs. [TBI and [T71 plot the resulting fit curves for each en- 
semble and the chiral/continuum extrapolated curve with 
its error band for fo(q 2 ) versus E\ (we show separately 
the coarse and fine ensembles in order to avoid too much 
clutter). The fit is excellent and has y^/dof = 0.44. In 
Table IVT1 we summarize fit results for * Df., k = 0,1,2 
coming from the simultaneous fit both for individual en- 
sembles and in the physical limit. These arc plotted in 
Figs. M - H3 In Fig. [UJ we plot f (q 2 = 0) for the 
five ensembles and in the physical limit. One sees that 
within errors this quantity shows little light quark mass 
dependence and a ~ 1.3% lattice spacing dependence. 

We call the chiral/continuum extrapolation based on 
the ansatz (J33J) - (EH) and shown in Figs. [Tol[T71 and[2T1the 
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FIG. 17: Chiral/continuum extrapolation of foil 2 ) versus 
E\ from the modified z-expansion ansatz. The data points 
are fine lattice points. Two individual curves and the extrap- 
olated band are from a fit to all five ensembles. 



FIG. 19: Same as Fig. [18] for ai * D x . 
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FIG. 20: Same as Fig. [18] for a 2 * D 2 . 



FIG. 18: The expansion parameter ao * Do versus the light 
quark mass from a simultaneous fit to all data. 



"simultaneous modified z-cxpansion extrapolation." We 
have tested the stability of this extrapolation by adding 
further terms to the ansatz and/or modifying some of the 
fit parameters and checking for changes in the physical 
limit fo(q 2 = 0). For example, we have, 

1. added xf terms, 

2. modified the lattice spacing dependent terms: 

(a) drop (am c ) 4 and (aEx) 4 terms 

(b) drop (am c ) 4 term 

(c) drop (aEx) 4 term 

(d) add (am c ) 6 term 



(e) add up to (am c ) 10 terms 

(f) add {aE K f term, 

3. used P{q 2 ) = 1 for the pole term, 

4. used Ek from dispersion relations, 

5. used an overall factor (1 + /(^SAI^ + SAIk)) outside 
the ^2 k to incorporate sea quark effects, rather than 
include them in the ZVs, 

6. used an overall factor (1 + d(am c ) 2 + e(am c ) 4 ) out- 
side the ^2 k to estimate the am c errors, rather than 
include them in the -CVs, 

7. replaced 6M 2 /K -> (M^ Tad ) 2 /(4^) 2 , 

8. used a simpler £)j = 1 + d\xi + d z 2 x s , 
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FIG. 21: fo at q 2 = for the five ensembles and in the 
physical limit. 
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FIG. 22: Tests of the "simultaneous modified z-expansion 
extrapolation" . The red horizontal line is the central value of 
the fit shown in Figs. I16l and ll7l and the orange lines indicate 
the error. The numbers under the data points correspond to 
the "test numbers" given in the text. 



9. used an even simpler Di = 1 + d\xi. 



Fig. [22] summarizes the results of these tests. One sees 
that the standard z-expansion extrapolation result is very 
robust. The second item of the tests checks that we es- 
timate the lattice spacing extrapolation error correctly. 
Until (am c ) 4 and (aEpc) 4 terms are included, the error is 
increasing, however after including the fourth powers the 
error is stabilized. This also shows that the am c error of 
the HISQ action is under control in our simulations. 



VII. RESULTS IN THE PHYSICAL LIMIT: /+(0), 
|V CS | AND UNITARITY TESTS 

This section summarizes the main results of this pa- 
per. We present our Standard Model prediction for the 
D -> K,lv decay form factor at q 2 = 0, /+(0) = /o(0), 
determine the CKM matrix element \V CS \ using input 
from BaBar and CLEO-c and carry out unitarity tests. 

/+(Q) - /o(Q) 

We have seen that the simultaneous modified z-expansion 
extrapolation method gives very stable results. It gives 
/+(0) = 0.748 ± 0.019 in the physical limit for D° -> 

K~lu, and /+(0) = 0.746 ± 0.019 for D+ -> K°lv. We 
take an average over these two channels and our final 
result in the physical limit becomes, 



f+~* K (fy = 0.747 ±0.011 ±0.015. 



(39) 



The first error comes from statistics and the second error 
represents systematic errors. Table IVIII summarizes the 
error budget. One sees that the largest contributions to 
the total error come from statistics followed by (am c ) 
and (clEk) extrapolation errors. 

In order to calculate the form factor, we have to put in 
meson masses from experiment and also from our lattice 
simulations. For example, we need experimental D, K, 
and 7T meson masses to get the form factor at the physical 
limit, and Ek, D, and K meson masses from the lattice 
calculations are used to fit at non-zero lattice spacing. In 
Table IVlTl "Input meson mass" refers to errors induced 
from these input meson masses. In the fit ansatz, eq. [Ml 
there are light quark (c\ and c|), strange quark (c\), and 
sea quark dependent terms (fi). Each systematic error 
due to these terms is shown on the fourth to sixth line 
in the table. Lattice spacing dependence errors are esti- 
mated separately for (am c ) n and [aEx) 3 type contribu- 
tions. 

In the fit ansatz, xilog(xi) is the most infrared sensi- 
tive term. We calculate the pion-tadpole loop integral 
both at finite volume and at infinite volume and com- 
pare these to estimate the finite volume effects. For the 
charm quark mass tuning error, wc calculate the form 
factor with a different charm quark mass, am c = 0.629, 
on the C3 ensemble, and compare with the result with the 
tuned am c = 0.6235 of Table HT1 All but the last two en- 
tries in Table [VlT] (finite volume and charm mass tuning) 
were calculated using methods introduced in reference 
[23l ] . The total error coming out of the chiral/continuum 
extrapolation can be decomposed into individual contri- 
butions, a 2 = J2i c i a ii where the sum goes over the 
first 8 entries in Table IVlTl Details are described in Ap- 
pendix B. 

One might worry about other potential systematic er- 
rors, not listed in Table IVlTl such as those due to missing 
sea charm quarks or electromagnetism/isospin breaking. 
The separate numbers given above eq. ([3"9"]> for D° K~ 

and D + — > K° form factors take into account just the dif- 
ferences in masses of the charged versus neutral mesons. 
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TABLE VII: Total error budget. 



Type 


Error 


O liCLlilo LlOcLl 


1 5 % 


IjcLL LlL-fc bCcLlc 1 / 1 eLIlCl / 1 / Ct J 


u.z /o 


Input meson niEtss 


u.i /o 


Light cnictrk dependence 


U.D /o 


?lfrflnp"p miflrV HptiptiHptipp 


0.7 % 


Sea quark dependence 


0.4 % 


am c extrapolation 


1.4 % 


aEpc extrapolation 


1.0 % 


Finite volume 


0.01 % 


Charm quark tuning 


0.05 % 


Total 


2.5 % 



2.5 % error 



"0- 



HPQCD (2010) 

Other Iheory 
Fermilab/MILC (2005) 

Sum Rules (2009) 

Experiment + CKM Unitarity 
CLEO-c (2009) 
BaBar (2007+update) 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

f + (q 2 =0) 



This "kinematic" effect is seen to be less than ~ 0.3%. 
It is much harder to assess the true dynamical electro- 
magnetic effects. However no statistically significant dif- 
ferences have been observed experimentally [f| and we 
will ignore further electromagnetic/isospin breaking ef- 
fects. Similarly we will assume that errors due to missing 
sea charm quarks are small enough so that they do not 
change the 2.5% total error when added in quadrature. 
This has been true in the case of several quantities where 
it was possible and appropriate to apply perturbative es- 
timates of dynamical charm quark effects po| . 

The total error for /+(0) is estimated here to be 2.5%. 
This is a factor of four times smaller than in the previ- 
ous lattice calculation of Ref. [(|. This was achievable 
because of applying several new methods and techniques 
that were described in the text. We employ the HISQ 
action for both charm and light quark actions and a 
scalar current rather than the traditional vector current. 
Because of these new methods, we obtain results with 
smaller discretization errors and no operator matching. 
We also developed the modified z-cxpansion extrapola- 
tion method, which is crucial to decrease errors due to 
the discretization, chiral / continuum extrapolation and 
parameterization of the form factor. In order to decrease 
statistical errors, we apply random-wall sources and per- 
form simultaneous fits with multiple correlators and X"s. 
If we compare with the error budget of Ref. Q , then we 
see the statistical errors reduced from 3% to 1.5% and 
the extrapolation and parameterization errors from 3% 
to 1.5% as well. The biggest improvement is in the dis- 
cretization errors. The total discretization errors have 
now been reduced from 9% to 2%. We note that the 
concept of the discretization errors is different in Ref. [f| 
compared to here. In Ref. @, they estimate the dis- 
cretization errors by power counting, since they calculate 
at only one lattice spacing. However we actually perform 
continuum extrapolations with correction terms for the 
discretization effects. As a result, we do not have dis- 



FIG. 23: Comparisons of /o(g 2 = 0) with other calculations 
and experiments. 

cretization errors per se, but instead extrapolation errors 
due to higher order correction terms. 

In their papers both BaBar Q and CLEO-c @ have 
converted their measurements of /+(0) * | V^ s | into re- 
sults for /+(0) using values for |V CS | fixed by CKM 
unitarity. For this CLEO-c uses the 2008 PDG CKM 
unitarity value of \V CS \ = 0.97334(23) [H| and obtains 
f^ K (0) = 0.739(9) and BaBar uses \V CS \ = 0.9729(3) 
leading to /+(0) = 0.737(10). In Fig. [23] we plot the 
new HPQCD result of this article, cq. ([39]) . together with 
earlier theory results from the lattice |6j and from a re- 
cent sum rules calculation [25[ and with the BaBar and 
CLEO-c numbers. One sees the very welcome reduction 
in theory errors which are now small enough so that the 
agreement between theory and experiment in Fig. [23] al- 
ready provides a nontrivial indirect test of CKM uni- 
tarity. We can, however, do better and carry out more 
direct tests of unitarity by determining |V CS | without the 
assumption of unitarity. 

Direct Determination of |V CS | 

As experimental input we take /+(0) * | V^ s | = 0.719(8) 
from CLEO-c @ and /+(0) * \V CS \ = 0.717(10) from 
BaBar [||. For the latter we have multiplied BaBar's 
quoted /+(0) with their quoted CKM unitarity value for 
\V CS \. Averaging between the two experiments we use 
/+(0) * \ V ca \ = 0.718(8) together with eq.(J3H]) to extract 
\V CS \. One finds, 

\V CS \ = 0.961 ±0.011 ±0.024, (40) 

in good agreement (as expected from Fig. |2"3"]) with the 
CKM unitarity value of 0.97345(16) 0. The first error in 
([40]) is from experiment and the second from the lattice 
calculation of this article. This is a very precise direct 
determination of \V CS \, made possible by the many ad- 
vances in lattice QCD that are described in this article 
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HPQCD (2010) 
l\g = 0.961(11)(24.) 

PDG 2010: direct estimation 
Semi-leptonic decay: IV J = 0.98 (1)(10) I — 
Leptonic decay: IV I = 1.030 (38) 
Average: IVJ = 1.023 (36) 

PDG 2010: CKM Unitarity 
IV 1 = 0.97345 (16) 
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FIG. 24: Comparisons of our new \V CS \ with values in the FIG. 25: Unitarity checks for the second row and second 
PDG 0. column of the CKM matrix. 



together with the tremendous progress in recent experi- 
mental studies of D semileptonic decays [Hill- In Fig. [24] 
we plot several previous direct determinations of \V CS \ 
from the 2010 PDG @ together with gDJ and the CKM 
unitarity value. 

In a companion paper [2(| where we update HPQCD's 
D s meson decay constant /d s , we also determine \V CS \ 
from D s — > tv and D s — > jiv leptonic decays. One finds 
I V cs I leptonic = 1.010(22) which is consistent with eq. (|40)) 
at the 1.4 a level. 

Further Unitarity Tests 

Using the new value of \V CS \, eq. (|4TJ]) . and the current 
PDG values \V cd \ = 0.230(11) and |V cf> | = 0.0406(13) one 
finds, 

\V cd \ 2 + \V CS \ 2 + \V cb \ 2 = 0.978(50) (41) 

for the 2nd row. And similarly for the 2nd column, with 
\V US \ = 0.2252(9) and \V U = 0.0387(21) one gets, 

\V US \ 2 + \V CS \ 2 + \V ts \ 2 = 0.976(50) (42) 

These 2nd row and 2nd column unitarity test results are 
shown in Fig. [55] together with the PDG numbers men- 
tioned already in the Introduction. Again one sees the 
improvement coming from the reduction in the uncer- 
tainty in \V CS \- 



These determinations serve mainly as consistency checks 
on our quark mass tunings and on our fitting and chi- 
ral/continuum extrapolation methods. More extensive 
studies of fo 3 involving five lattice spacings are reported 
in [2(|. Here we present the first results for f D us- 
ing HISQ charm and light quarks that employs the new 
HPQCD n scale [l9[ (the scale we use throughout in this 
article) . 

We use continuum partially quenched ChPT formulas 
augmented by discretization terms to extrapolate to the 
chiral/continuum limit. Error budgets are close to those 
in ref fiol ] with, however, a decrease in the n uncertainty. 
For f n and fx, we find, 

U = 132.3 ±1.6 MeV (43) 

Sk = 157.9 ± 1.5 MeV (44) 

One sees that agreement with experiment is good and 
within 1 - 1.5o- (or 1.5%). 

For fjj we have results from the simultaneous fit with 
three-point correlators and also from just the two-point 
correlator fits. As we noted in section V, the former 
has smaller errors however they are consistent with each 
other. 

jSimFit = 2Q6 3 ± 4 3 Me y ( 45 ) 



VIII. FURTHER RESULTS FROM TWO-POINT f P t = 2 n i ± 5 7 MeV (46) 

CORRELATORS D { ' 

Both values show good agreement with experiment as 

In this section we summarize physics results extracted can be seen in Fig. [26] 

from two-point correlators that emerged as part of our For fjj , we find, 
analysis of D meson semileptonic decays. This includes 

determinations of decay constants, f n , Jk, Id and Jd s - Id, = 250.2 ± 3.6 MeV, (47) 
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FIG. 26: Comparisons of fo with experiment, simultaneous 
fit, and two-point correlator fit. 

which agrees with (but is less precise than) our recent 
update in reference [20j. fn, = 248.0 ± 2.5 MeV. One 
sees that with experimental values having come down in 
recent years and with the increase in the HPQCD value, 
there is no longer any discrepancy (beyond 1.6 a) between 
theory and experiment. The current HFAG [26| number 
is f Ds = 257.3 ±5.3 MeV. 
Finally we present the ratio, 

+ - — = 2.986 ± 0.087 GeV' 1 . (48) 

JD S 

This quantity can also be obtained from experimental 
measurements of D — > K 1 Iv scmilcptonic and D s leptonic 
decay branching fractions, and has the virtue that |V CS | 
drops out in the ratio. We compare (|48|) with experiment 
in Fig. [57] and good agreement is found. 

IX. SUMMARY AND FUTURE OUTLOOK 

We have completed the first study of D — > K semilcp- 
tonic decays using the HISQ action for the valence charm, 
strange and light quarks. The most important result of 
this article is given in eq. Q39p and provides the form fac- 
tor f+~^ K (q 2 ) at q 2 = 0. We were able to determine 
this quantity with a 2.5% total error which represents a 
four fold improvement in precision over earlier theory re- 
sults. This is shown in Fig. 1231 We then combined our 
form factor result with recent measurements of D — > K 
scmilcptonic decays by the BaBar and CLEO-c collabo- 
rations to extract a very accurate direct determination of 
the CKM matrix element | V cs \ . This is given in eq. P0|) 
and comparisons with previous determinations shown in 
Fig. [Ml The new value for |V CS | is consistent with the 
PDG value based on CKM unitarity. We carried out 
direct tests of 2nd row and 2nd column unitarity with 



FIG. 27: Comparisons of f^ K (0)/f Ds . We estimate the 
experiment value from a simple calculation of combining fr> s 
from HFAG [H and /+^ A '(0) from BaBar [| and CLEO- 
c H. This is an appropriate thing to do since both quantities 
were derived from experimental measurements using the uni- 
tarity value of | V cs | ■ 



results given in equations (|41[) and (|4"2"j) and depicted in 
Fig. [23 Although still far from the accuracy achieved 
in 1st row unitarity tests, the reduction in errors on \V CS \ 
has made 2nd row and column unitarity tests much more 
relevant and interesting than in the past. In section VIII 
we give a lattice QCD value for f+^ K {0)/fD B which is 
consistent with experiment. Within the current theory 
and experimental errors, this provides a highly nontriv- 
ial consistency check on how we treat D scmilcptonic and 
leptonic decays on the lattice and more generally in the 
Standard Model. 

The calculations of this article can be improved upon 
and extended in several ways. The largest errors in 
Table I VIII from statistics and continuum extrapolations 
can be reduced straightforwardly by working with more 
gauge configurations and time sources and at more values 
of the lattice spacing. An obvious extension of the cur- 
rent study is to investigate D — > ir scmilcptonic decays 
and determine \V c d\. Work on this project has already 
begun. In a future project we also plan to calculate the 
vector current hadronic matrix elements {K{-k)\V^\D) . 
This will provide f+(q 2 ) as a function of q 2 . As men- 
tioned in section II we will be carrying out nonperturba- 
tive matching of the vector current based on PCVC for 
this project. 
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TABLE VIII: A sample set of the priors and prior widths 
for ensemble Fl with p= (0,0,0) and p = (1, 1,0). We have 
tested with various priors and prior widths, and the fit re- 
sults are not sensitive to reasonable variations. Note that 
i = 1,2, 3, ... and j, k = 0, 1, 2, .... 





prior 


width 


prior 


width 




p= (0,0,0) 




P= (i,i,o) 




A 


n m 


D 1 


U.U-L 


D 1 


TD 

-Djfc 


U.U1 


U.l 


U.U1 


O 1 

U.l 


Gjfc 


U.U1 


U.l 


n ai 
U.U1 


U.l 




0.01 


0.1 


0.01 


0.1 


fP 


SI ^ 

U.O l-O 

U.O 


+0.815 

-0.408 
+0.6 
-0.3 


U.O 

n a 
U.4 


+0.8 

-0.4 
+0.4 

-0.2 


bE 


0.12 


1.0 


0.015 


0.3 


b D 


0.1 


1.0 


0.03 


0.3 


rplD t?ID 


i n 
0.4 


+ 1.0 
-0.5 
+0.4 
-0.2 


l.U 

0.4 


+ 1.0 

-0.5 
+0.4 
-0.2 




0.01 


0.1 


0.0028 


0.1 


d? 


0.01 


0.1 


0.006 


0.1 


Eg 


0.23 
0.5 


+0.23 
-0.12 
+0.5 
-0.25 


0.39 
0.4 


+0.39 
-0.2 
+0.4 

-0.2 


bf 


0.15 


1.0 


0.01 


0.1 


Eq 

&i — 


0.4 
0.5 


+ 0.4 
-0.2 
+0.5 
-0.25 


0.53 
0.4 


+0.53 
-0.27 
+0.4 
-0.2 




0.01 


0.1 


0.001 


0.01 
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Appendix A: Priors and Prior Widths for Section V 

We list a sample set of priors and prior widths in Ta- 
ble IVIIII that have been used for the simultaneous fits 
to Cj^ nt , C^ nt , and three of the three-point correlators 
with T = 19, 20, and 23 for ensemble Fl. The fit ansatz 
and results have been presented in Section V. 



Appendix B: Bayesian Fits in Section VI 



From ([33)1 - (|38[) one sees that the chiral/continuum 
extrapolation ansatz for fo(q 2 ) starts out with 23 basic 
fit parameters c„, n = 1,2,. ..,23. 



C n ■ an, Ol, Cl2, Ol, o 2 , c\, C 2 , C3, di, Ci, fi. 



(Bl) 



where i = 0, 1, and 2. In a Bayesian fit each of these 
fit parameters will have its own prior c„ and prior width 
a n , which we call "Group I." Our choices for these priors 



will be discussed below. Our fit ansatz for /n includes 
in addition to the fit parameters c„ also many input pa- 
rameters such as Md, Ek, J"i etc. all of which have some 
uncertainty associated with them. Ref. [23| describes a 
method to include effects coming from these types of un- 
certainties into the final error in the extrapolated value 
for /n. What one does is convert all these input param- 
eters into additional fit parameters with priors and prior 
widths given by their known central values and errors Us- 
ing this approach we have changed 61 input parameters 
into new fit parameters pj, j = 1,2, ....,61, which we call 
"Group II." 

p 3 : {^)\aAP Dl aE l K {p),aMl,aMl, 
{aMp tad )\ (aMZ sqtad )\M l D , , 
n , M^ ys , Mf ys , Mf ys , Mf ys , M p *r (B2) 

where i = 1,2, .., 5 goes over the 5 ensembles. The use of 
Group II parameters is a very efficient way to include er- 
rors coming from uncertainties in input parameters into 
our final total error. An alternative approach would re- 
quire visiting each input parameter in turn, redoing the 
chiral/continuum extrapolation and coming up with an 
estimate for the systematic error coming from this input 
parameter. In our approach no additional systematic er- 
rors for these input parameters are called for and the 
effects from their uncertainties are included in the ex- 
trapolation error. For instance, this is very helpful to 
include the error from E\. In Fig. [TBI and [171 there are 
errors for E\ on the lattice data as well as the extrapo- 
lated results, and estimating these errors is not a trivial 
task. However, introducing the Group II parameters in- 
corporates all E 2 K errors as part of the final vertical error. 

In Bayesian fits one minimizes the augmented 
chisquared, 



A a 



Xtr aditional 



^Group I 



XGroup II 



2 1 2 1 2 

Xtr aditional roup 1 "+" Xg roup II 
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(Jo - foiansatz)) 2 

Ko) 2 



E 

i=l 
2:5 

E 

n=l " n 

L> a 2 

3 = 1 3 



(B3) 
(B4) 

(B5) 

(B6) 



When carrying out the chiral/continuum extrapola- 
tions we have expressed all dimensionful quantities in 
units of GcV. We give the set of priors and prior widths 
for the Group I parameters c n in Table IIXI and for the 
Group II parameters pj in Table [X] and IXII 

Setting priors and prior widths is straightforward. For 
Group I, the quark mass terms, such as xi and x s , are nor- 
malized by the scale, A = AttJ^. Therefore, it is natural 
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TABLE IX: Priors and prior width of the Group I parameters 
for the simultaneous modified ^-expansion extrapolation fit 



Group I 


prior 


width 


fit result fit error 


do 





1 


0.09766 


0.0029 


ai 





1 


0.08999 


0.02 


a 2 





1 


-0.07044 


0.11 


bi 





0.3 


0.03775 


0.13 


b 2 





0.3 


0.07179 


0.17 


c? 





1 


-0 52596 


0.31 


ci 





1 


-0.19051 


0.82 


ci 





1 


0.02877 


1 







1 


-0 09919 


0.98 


c\ 





1 


0.00827 


1 


4 





1 


0.00044 


1 


4 





1 


-0.02897 


0.24 


4 





1 


0.32804 


0.66 


4 





1 


-0.03116 


1 


do 





0.3 


0.00966 


0.11 


di 





0.3 


0.01769 


0.29 


d 2 





0.3 


0.00541 


0.3 


eo 





0.2 


0.01554 


0.19 


ei 





0.2 


0.00447 


0.2 


ei 





0.2 


0.00131 


0.2 


fo 





0.3 


-0.10860 


0.28 


fi 





0.3 


-0.00474 


0.3 


h 





0.3 


0.00105 


0.3 



that we expect that the parameters vary between —1 to 
1. However, we know from other lattice calculations with 
the same gauge configurations and our lattice data that 
the sea quark mass contribution is smaller than that of 
valence quark. Thus, we take the priors and prior widths 
for the sea quark mass terms as fi = ± 0.3. The lead- 
ing heavy quark error is proportional to 0(a a (am c ) 2 ) 
so we use dj = ± 0.3 for the (am c ) 2 terms. For the 
purposes of setting priors, we conservatively do not in- 
clude a factor of v 2 /c 2 here. On the other hand for the 
(am c ) 4 terms we do take the expected factor of v 2 /c 2 
into account and choose &i — ± 0.2. Similarly we use 
b 3 ■ = ± 0.3 for the (aE K ) 2 and (aE K ) 4 terms. This 
reflects a factor of a s for the (aEK) 2 terms and the fact 
that higher powers of (ap M ) typically come with smaller 
numerical factors relative to lower powers (such as in an 
expansion of ^sinh(aE)). For Group II, we use lattice 
results and experiments that we described in the text for 
the priors and prior widths. 

As we stated above, all sources of systematic errors 
are already included in the fit ansatz, except the finite 
volume and charm quark mass tuning errors. We can 
consider the total error squared, a 2 , as a linear combina- 



TABLE X: Priors and prior width of the Group II parame- 
ters for the simultaneous modified z-expansion extrapolation 
fit. Parameters with five rows correspond to that on the five 
ensembles, CI, C2, C3, FI, and F2. 



Group II 


prior width 


fit result fit error 


r\ 

M phya 
M vhya 
M Phys 

M phys 


0.3133 0.0023 
0.6858 0.004 
0.1373 0.0023 
1.8645 0.0004 
0.4937 0.000016 
2.3173 0.0006 


0.313285 0.0023 
0.685799 0.004 
0.1373 0.0023 
1.8645 0.0004 
0.4937 0.000016 
2.3173 0.0006 


ri/a 


2.647 0.003 
2.618 0.003 
2.644 0.003 
3.699 0.003 
3.712 0.004 


2.64677 0.003 
2.61818 0.003 
2.64388 0.003 
3.69905 0.003 
3.71213 0.0039 


aM D 


1.13927 0.00066 
1.15947 0.00076 
1.16179 0.0005 
0.814062 0.00035 
0.819663 0.00026 


1.13925 0.00066 
1.15949 0.00076 
1.16179 0.0005 
0.814066 0.00035 
0.819663 0.00026 


aE K 
P= (0,0,0) 


0.312174 0.0002 
0.32851 0.00048 
0.357205 0.00022 
0.228546 0.00017 
0.24596 0.00014 


0.312159 0.0002 
0.328529 0.00048 
0.357228 0.00022 
0.228572 0.00017 
0.245944 0.00014 


aE K 
p= (1,0,0) 


0.408141 0.00066 
0.453061 0.0016 
0.475036 0.00087 
0.320348 0.00069 
0.333981 0.00037 


0.408059 0.00065 
0.453098 0.0015 
0.474849 0.00086 
0.320396 0.00067 
0.333932 0.00036 


aE K 

P= (1,1,0) 


0.483702 0.00094 
0.552504 0.0017 
0.572004 0.001 
0.391932 0.00087 
0.401385 0.00065 


0.483798 0.00093 
0.552527 0.0016 
0.572016 0.001 
0.391883 0.00085 
0.401619 0.00064 


aEK 
P= (1,1,1) 


0.546941 0.002 
0.637288 0.0032 
0.652424 0.0022 
0.45593 0.0015 
0.460897 0.0011 


0.547535 0.002 
0.637096 0.0032 
0.652314 0.0022 
0.455311 0.0015 
0.460966 0.0011 


aM, l3 


0.411128 0.00018 
0.414346 0.00022 
0.411848 0.00022 
0.294159 0.00012 
0.293114 0.00018 


0.41113 0.00018 
0.414345 0.00022 
0.411847 0.00022 
0.294158 0.00012 
0.293117 0.00018 




0.159893 0.00017 
0.210815 0.00023 
0.293124 0.00023 
0.134449 0.00015 
0.187346 0.00013 


0.159895 0.00017 
0.210813 0.00023 
0.293119 0.00023 
0.134445 0.00015 
0.18735 0.00013 
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TABLE XI: (Continued) Priors and prior width of the Group 
II parameters for the simultaneous modified z-expansion ex- 
trapolation fit. Parameters with five rows correspond to that 
on the five ensembles, CI, C2, C3, Fl, and F2. We quote 
M* sqTad and M? sqTad from ref. M- 



Group II 


prior width 


fit result fit error 




0.3653 0.00029 
0.38331 0.00024 
0.40984 0.00021 
0.25318 0.00019 
0.27217 0.00021 


0.365304 0.00029 
0.383309 0.00024 
0.409839 0.00021 
0.253177 0.00019 
0.272174 0.00021 


aM aaqtad 


0.15971 0.0002 
0.22447 0.00017 
0.31125 0.00016 
0.14789 0.00018 
0.20635 0.00018 


0.15971 0.0002 
0.22447 0.00017 
0.31125 0.00016 
0.147889 0.00018 
0.206351 0.00018 




2.32897 0.005 
2.32899 0.005 
2.33053 0.005 
2.32293 0.005 
2.32113 0.005 


2.33033 0.0047 
2.32818 0.0048 
2.33028 0.0048 
2.32179 0.0045 
2.32197 0.0045 



tion of each source [2S 



2.-! 

E 

n=l 



61 



Cr 



£3 



(B7) 



where the first term is for Group I (c„), and the sec- 
ond term is for Group II (pj). We actually calculate the 
contributions from each source, C Cn a 2 n and C Pj <r 2 . us- 
ing the method presented in [23j j , and they add up to the 
total error a 2 correctly. We group together appropriate 
parameters, and list them in Table IV*lT1 



Appendix C: Chiral and Continuum Extrapolations 
based on Chiral Perturbation Theory 

In this Appendix we carry out further consistency tests 
of the chiral/continuum extrapolation of section VI by 
working with a completely independent fit ansatz. We 
will use the partially quenched chiral perturbation the- 
ory (PQChPT) formulas developed in Refs.flol. ITm] aug- 
mented by terms parameterizing discretization effects 
and Ek dependence. 

Heavy meson ChPT formulas arc organized through 
form factors f\\ and f± in terms of which /o(<? 2 ) is given 
by, 



fo(q 2 ) 



V2M L 



Ml 



Ml 



[(M D - M K )f {l + (E 



M 2 K )f±] 



1.1 



; 0.9 



0.8 



0.7 



chiral/conti. extrapolation 
C3 ensemble 
C2 ensemble 
CI ensemble 

f„atq 2 = 




0.2 



0.4 



0.6 



2.4% error 

_J I I I I 

0.8 1 1.2 



E K [GeV 



FIG. 28: Chiral/continuum extrapolation of /o(<7 2 ) versus 
E 2 K from the ChPT ansatz. The data points are coarse lattice 
points. Three individual curves and the extrapolated band are 
from a fit to all five ensembles. 



We follow very closely the approach and notation of 
Ref. [13], however with all the taste breaking effects 
turned off. /n and f± are parameterized as, 



./'ll 



K 



1 + <5/|| + cjmi + c^m S ' + 

4 ea (2m a + m s ) + h\\ {E K ) (1 + c (am c ) 2 
+ci(am c ) 4 ) (C2) 
[l + 5f±+clmi + 



fjn 



(CI) 



U(E k +A*+D) 
(cf + cf- - c\,)m s , + cj ea (2m u + m s ) + 
h±(E K )] (1 + c (am c ) 2 + ci(am c ) 4 ). (C3) 

The chiral logs are contained in 6f\\, Sf± and D. We 
give their explicit expressions in Appendix D. mi and 
m s i arc valence and m u and m s the sea quark masses. 
<?7r is the DD*tt coupling and A* the D* — D mass 
splitting. h^(Ex) and h±(Ej() are unknown functions 
of Ek- We will use polynomial expansions for them. 

h\i±{E K ) = c\^E K + ctt^K + Tnc first tw0 tcrms 
are motivated by ChPT jl7l]. however we view h\\ t j_(EK) 
as potentially parameterizing Jo more generally, even be- 
yond the regime of small Ex where ChPT is assumed 
valid. Our data is fit very well (x 2 /do/=0.48), however, 
all the way to Ek ~ lGeV keeping just terms through 
0{E K ) in h\\ and hj_. Figs . l28l and [29l show results from 
a simultaneous fit to all our data points using the ChPT 
ansatz. These should be compared to Figs. [TH] and [T7] 
from the modified z-expansion ansatz. 

In Fig. [30] we compare fo(q 2 ) in the physical limit com- 
ing from the z-expansion extrapolation of section VI and 
the ChPT extrapolation of this Appendix over the entire 
physical q range. And in Fig. [21] we compare results at 
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FIG. 29: Chiral/continuum extrapolation of /o(<? 2 ) versus 
E\ from the ChPT ansatz. The data points are fine lattice 
points. Two individual curves and the extrapolated band are 
from a fit to all five ensembles. 



FIG. 31: Comparisons of /o at q 2 = for the five ensembles 
and in the physical limit from the z-expansion and the ChPT 
extrapolations. 



Appendix D: Partially Quenched ChPT Chiral Logs 



1.1 



; 0.9 



0.8- 



0.7 



— Simultaneous ChPT fit 
modified z-expansion fit 

O f at q = (ChPT) 

9 
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0.4 



0.6 
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FIG. 30: Comparisons of fo(q 2 ) in the physical limit from 
the z-expansion and the ChPT extrapolations. 



q 2 = for each ensemble and in the physical limit. One 
sees that the two extrapolations are nicely consistent with 
each other. We believe the consistency check of this Ap- 
pendix has been very useful. It provides further support 
for the results of the z-expansion extrapolation and in- 
dicates that errors there were not underestimated. Note 
that the ChPT ansatz includes all the complicated chi- 
ral logs of Appendix D, whereas the z-expansion ansatz 
makes do with just a simple xilogxi term. 



In this appendix we summarize partially quenched 
ChPT (PQChPT) expressions for the chiral logarithm 
terms 6f\\, Sf± and D that we employ in cqs. (|C2[) and 
(jC3l) . The formulas for PQChPT in continuum QCD for 
heavy-to-light semileptonic decays were first developed 
in Ref.[l6| for degenerate sea quarks. Ref. 17 1 general- 
ized these results to nondegenerate 1+1+1 sea quarks 
and also to Staggered ChPT. We have started from the 
1+1+1 continuum PQChPT expressions given in Ref. 
[I?} for individual diagrams and for the D and kaon wave 
function renormalizations to obtain the full hi the 

2+1 PQChPT case. For the convenience of the reader 
we give these reconstructed expressions below. We use 
the same notation as in Ref. [T7| . u x" and "y" stand for 
the light valence quarks in the daughter meson (for the 
kaon x = I and y = s') and "it" and "s" denote sea 
light and strange quarks. Furthermore m a b is the mass 
of the pseudoscalar meson with quark content a and b 



and 



h(m 2 uu + 2m 2 ss ). 



(47r/) 2 r5/, 



D^K 



•3^ 



h{m xu ) + -h (m xa ) 



\(m V u) + 2 J l(«V) 
[2I 2 (m yu ) + h(m ya )] 

XX ) (h(m xx ) + h{rn xx )) + 
Rf 2] {m yy ) (h(m yy ) + h{m yv )) + 
Rf 2] (m n ) (h(m v ) + I 2 (m v )) 

^■[DR^m^h) 



1 ^DR^{m yy ;h) 
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+ -\DR^\m yy] h) 
(Anf) 2 6f^ K = 



h(m yu ) + -h{m ys ) 



-3^ 



(Dl) 



h{m xu ) + -h(m xs ) 



Rf 2 \m) 
R^(m;l) 



( m lu - 


■m 2 )(m 2 s - 


TO 2 ) 


( m uu 


~ m )\ m ss 


TO 2 ) ' 
— TO 2 ) 




(m 2 — to 2 ) 




(m 2 uu 


-m 2 )(m 2 sa 


-TO 2 ) 


(to 2 — TO 2 ) 



(D10) 



T(m) 



R§' 2] (m xx ) K 1 {m xx )+ 
Rf 2 \m yy )K x { 



l(m v ), 

(Dll) 



and 



DR^ 2 \m yy -h) 



6 



DR?> 2 \m xx -h) 



DR^ 2 \m;l) = ^-R^ 2 \ m; l). 
om 



(D12) 



(D2) 



{AvffD 



2 r~,D^K 



-3gl( v -p) x 



[2K x {m yu ) + K^mys)] + - \DR^(m yy ;Kx] 

(D3) 

In the D meson restframe v ■ p = Ek ■ Furthermore one 
has, 



Ii(m) = m 2 log 



A 2 



(D4) 



J a (m) - -2(v ■ vflog^ - i(v ■ P ) 2 F f ) + 2(v ■ p) 2 



v ■ p 



with 

F(x) = 



(D5) 



VI - x 2 tanh- 1 ^! - x 2 ) < x < 1 



Vz 2 - 1 tan 1 (\ / x 2 - 1) x > 1 



(D6) 



Xi(m) = 



+ 3^ 



J 2 + 3 ^ ' P) 



2 21 

TO 



TO \ 10. , 2 4 2 2tt to 3 

* F( )-^(v-p) 2 + -m 2 -— , 

v -p) 9 3 3 v ■ p 



R [s,2] ,n = ( m L -™ 2 )(™L ~ W 2 ) 

X ( m yy — m 2 ){m 2 — to 2 ) 

#[3,2] / m ) = (m 2 M -?71 2 )(TO 2 s -TO 2 ) 
^ ( m xx — TO 2 ) (to 2 — TO 2 ) ' 



(D7) 
(D8) 

(D9) 



We refer the reader to the original literature [lg, [l7[ for 
further details. Here, for completeness, we give partially 
quenched formulas for the chiral logarithms in D — > tt 
decays. They will be used shortly in our own studies of 
D — > 7r, Iv decays. Some care is required in taking the 
y ->• x limit of (fPT]) . (|D2j) and (|D3|) . however in the end 
expressions are simpler for D — > tt than for D — > K . 



(4nf) 2 Sf^ 
(1-35 2 ) 



h(m xu ) + -zh{m xs ) 



[21 2 {m xu ) + h{m xs )} 



1+3^ 
6 



2 r 



DR^ 2 \m XXl h 

(D13) 



"(1 + 3^) 



) + -h(m 



DRp' 2 \m xx ;Ki) 



l + 3 5 2 



Dfll 2 ' 2 ]^;/!) 

(D14) 



(47r/) 2 ^— = -3g 2 .(vp) x 
[2A'i(to. tu ) + K\ (m xs )] + l \DR^ 2 - 2 \m xx - K\ 



(D15) 
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